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galactic halos. Assuming that the predicted dark matter clumps are composed of annihilating neu- 
tralinos, we examine their detectability via Atmospheric Cherenkov Telescopes (ACTs). Depending 
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on their density profile, individual neutralino clumps should be observable via their 7-ray continuum 
and line emissions. We find that the continuum signal is the most promising signal for detecting 



a neutralino clump, being significantly stronger than the monochromatic signals. Limits from the 
2 ' line detectability can help lift degeneracies in the supersymmetric (SUSY) parameter space. We 

show that by combining the observations of different mass clumps, ACTs can explore most of the 
SUSY parameter space. ACTs can play a complementary role to accelerator and 7-ray satellite 

X' 

limits by exploring relatively large neutralino masses and less concentrated clumps. We develop 
a strategy for dark matter clump studies by future ACTs based on VERITAS specifications and 
encourage the development of techniques to identify primaries. This can reduce the background 
by an order of magnitude. 
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I. INTRODUCTION 



The presence of a dark matter component is inferred through gravitational interactions 
in galaxies and clusters of galaxies. Its contribution is estimated to be about 30% of the 
critical density of the Universe. This dark matter cannot be all baryonic. Constraints 
from primordial nucleosynthesis and cosmic background radiation measurements limit the 
baryonic content of the Universe to be at most about 4% of the critical density. The 
remaining ~ 26% of the critical density is believed to be composed of a, yet to be observed, 
Cold Dark Matter (CDM) particle. 

Among the several candidates proposed for the non-baryonic dark matter, the leading 
scenario involves Weakly Interacting Massive Particles (WIMPs). Weakly interacting relics 
from the early universe with masses from some tens of GeV to some TeV can naturally give 
rise to relic densities in the range of the observed dark matter density. WIMPs are also well 
motivated by theoretical extensions of the standard model of particle physics. In particular, 
the lightest supersymmetric particle (LSP) in supersymmetric extensions of the standard 
model can be a WIMP that is stable due to R-parity conservation. In most scenarios the 
LSP is neutral and is called the neutralino, x (f° r a review see, e.g., U). 

Neutralinos may be detected directly as they traverse the Earth, may be generated in 
future particle colliders, and may be detected indirectly by observing their annihilation 
products. Direct neutralino searches are now underway in a number of low background 
experiments with no consensus detection to date. Collider experiments have been able to 
place some constraints on the neutralino parameter space (see, e.g., 0]), but a large region of 
the parameter space is still unexplored. Indirect searches offer a viable and complementary 
alternative to direct searches and accelerator experiments. In particular, neutralino annihi- 
lation in the galactic halo can be observed via 7-rays, neutrinos, and synchrotron emission 
from the charged annihilation products [g § |, |, |7|, g, g 0, [llj Q . 



The rate of neutralino annihilation is proportional to the neutralino number density 
squared (~ n 2 ), thus the strongest signals are likely to come from either the galactic center || 
H> H) @> 0) HI) or the higher density nearby clumps of dark matter in our halo |8|, ^ [10], |TTJ. 
The dark matter density in the galactic center depends strongly on the formation history 
of the central black hole ||13|| , and its 7-ray emission due to neutralino annihilation may be 
difficult to distinguish from other baryonic 7-ray sources. As an alternative, the clumpy 
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nature of CDM halos provides a number of high density regions, other than the galactic 
center. 

Clumps of dark matter distributed throughout galaxy halos is a generic prediction of high- 
resolution CDM simulations (see, e.g., [13, |Tj| |Tj| |T7[]). This high degree of substructure is 



what has remained from the merging of small halos to form the larger halos (either cluster or 
galaxy halos). Simulations show that about 10% to 15% of the total mass of a given halo is 
in the form of smaller mass clumps. These clumps are ideal sites for neutralino annihilation 
and can be observed via their 7-ray signal || ||, |10| , and via their synchrotron emission [ ] I 



Here, we study in detail the 7-ray emission from the nearest clumps of dark matter via 
Atmospheric Cherenkov Telescopes (ACTs), such as VERITAS @, HESS @, MAGIC gO 



and CANGAROO-III |2T[. We explore the neutralino parameter space for the continuum, 
the 77, and the Z7 lines and calculate the flux of different mass clumps assuming different 
clump density profiles and VERITAS specifications. We find that the larger mass clumps 
can be easily detected by the next generation ACTs and that the detectability of the smaller 
mass clumps can help constrain the neutralino parameter space. Future ACT observations 
are only limited by their threshold energy and will be more effective for neutralino masses 
above ~ 50 to 100 GeV, a range of neutralino masses that nicely complements accelerator 
and satellite studies of the neutralino parameter space. 

In the next section, we discuss the structure and the mass spectrum of the dark matter 
clumps that we assume in our study. In §111, we present the neutralino parameter space 



for continuum 7-rays, the 77, and the Z7 line emission using DarkSUSY p8| . In §IV, we 
discuss the characteristics of ACTs and the relevant backgrounds. Our results and strategies 
for CDM clump detection are shown in §V. Finally, we conclude in §VL 

II. THE DARK MATTER CLUMPS 
A. The nearest clumps 



According to N-body simulations 0, 15), pL6|| , substructures of masses > 1O 7 M survive 



in significant numbers in galactic halos. Within the virial radius of a galactic halo about 
10%-15% of the mass is in the form of dark matter clumps. The mass spectrum for the 
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simulated dark matter clumps in galactic halos is of the form (e.g., |15, |16| ) 

dN dump a 

dM dump ^ M ^ ' Uj 

for clumps with mass M dump ^10 7 M & , the lowest resolvable mass, approximately. Clumps 
are more numerous the less massive they are, thus one would expect a large number of 
relatively light clumps. Simulations find a ~ 1.9 (e.g., ||T^[). We adopt the same value here 
and extrapolate to lower masses than those presently resolved by simulations. With this 
extrapolation we can estimate the numbers of clumps in specific mass bins, as well as the 
distances to the nearest of these clumps. 

In order to find the proportionality constant in Eq. ([!]), we normalize the mass function 
so that 500 satellites exist within the galactic halo with masses above 10 8 M Q , in agreement 



17 / -s ir 7, r \ rr\r\ \ lvl max lvl min) x , n \ 

N dump {M > M min ) = 500 x — 7^77 9 (2) 



with the predictions of N-body simulations for a halo like ours flip] . Thus we find for the 
predicted number, N dump , of clumps with mass M greater than M min , in solar masses, 

\M T max I -M ^miri) 1 

(M max /1O 8 M ) 

where M max is the maximum mass clump in the halo. Using our halo as the host within 
which the clumps are orbiting, we assume M max ~ l%Mh ost — 2 x 1O 1O M . 

The number, NM c[ump , of clumps predicted in several mass bins that lie in the mass range 
from 10 2 M Q up to 1O 8 M are tabulated in the first and second columns of Table |. We choose 
each bin labeled by M dump to contain clumps with mass in the range [0.5M dump , 5M dump ], 
and calculate N Mclump = N dump (M > 0.5M dump ) - N dump (M > 5M dump ). The numbers 
appear to be high enough, and even though they are consistent with simulations, the question 
whether they can be realistic is raised. Contrary to what was initially believed, namely that 
large numbers of substructures within a galactic halo might threaten the stability of the 
disk, nowadays we know that this is not the case (for a recent review, see fTTfl). It has 
been found, first, that the orbits of satellite substructures in present-day halos very rarely 
take them near the disk ||22|| , and second, that disk overheating and stability would become 
real problems only when interactions with numerous, large satellites (e.g., the LMC) take 
place |23], |24j| . Consequently, these numbers can be considered perfectly realistic. In the 
third column of Table | we give values of the mean inter-clump distance in a specific mass 
bin, assuming that the clumps are homogeneously distributed within the virial extend of the 

1 /3 

halo. Thus, in the third column we have the quantity d mean = R viria i/Nj^ ciump , with R V i r i a i 



TABLE I: Numbers of clumps, mean and nearest clump distances. 
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10 7 


6.6 x 10 3 16023 
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10 8 


8.2 x 10 2 


32013 


3874 



the virial radius of the galactic halo which we take to be equal to 300 kpc. However, clumps 
are unlikely to be homogeneously distributed in the host halo. It is believed that they are 
distributed in a way similar to the way the smooth component of the dark matter halo is 
distributed (e.g., p5| ). Thus, we estimate another measure for the distance, the distance to 
the nearest clump for each mass bin, as follows: first, we calculate the fraction, /, of the 
total dark matter of the galaxy which is in the form of clumps of mass M c i ump and take 
/ ~ N Mclump M dump /M host . For a local dark matter density ~ O.OlM /pc 3 and using the 
above fractions we find the local dark matter density which is in clumps of a specific mass. 
Going from the mass to the number density, we find the distance to the nearest clump, 
dnearest — 36kpc/NM cl (see also p6|, p7fl). As we can see from Table |, d nearest can be a 
factor of ten smaller than the mean distance as calculated in the homogeneous distribution 
case. We use d nearest below when estimating 7-ray fluxes from different clumps. 

One should keep in mind that since we only have one halo to probe, our galactic halo, 
the actual distances to the nearest clumps will have large fluctuations. For a number of 



realizations of the clumpy halo in 7-rays, see Jul . Instead of simulating the sky, we focus on 
the strategy for future ACT studies based on d near est- As will be shown, the results we derive 
can be easily extended to accommodate different clump distances from the ones we use. A 
deeper understanding of the clump space distribution awaits higher resolution simulations. 
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B. The structure of the dark matter clumps 



We model the clumps using three different mass density profiles, p(r), the Singular 
Isothermal Sphere (SIS), the Moore et al. [|HJ, and the Navarro, Frenk and White (NFW) |28j 
density profile. The SIS profile is described by 

P(r) = % (3) 

and represents a reasonable upper limit to the degree of central concentration of the clumps. 
At the center, the SIS density profile diverges as r~ 2 . The Moore et al. density profile is 
given by 

P( r ) = 1~, \i.5 F, . , i TOT (4) 
ir/r s ) [l + (r/r s ) \ 

where p Q and r s are the characteristic density and scale radius of the configuration, respec- 
tively. This profile is the one suggested by numerous N-body simulations and is used in our 
study as representative of the intermediate concentration cases, at least compared to the 
other two density profiles. It is divergent at the center where it behaves as r -L5 , whereas at 
large distances (r > r s ) it behaves as r -3 . The NFW profile |28| can be written as 



p{r) = P - ^ • ( 5 ) 

(r/r s ) (1 + r/r s ) 

This profile was the first proposed density profile for the dark matter halos produced in 
simulations. It is divergent at the center with an r _1 behavior, namely it is shallower, less 
centrally concentrated than both the SIS and the Moore et al. profiles. At larger radii, 
NFW behaves as r -3 , exactly the same way as the Moore et al. profile. Note that the 
two quantities, p Q and r s , in the NFW and Moore et al. profiles are not the same for the 
same clump (for a comparison, see P§fl). We choose to use the NFW profile as the least 
centrally concentrated density profile. The density profile derived from Taylor and Navarro 
(TN)[^, |3lf , which is even shallower at the center behaving as r~ 75 , is a possibility; it will 
give lower fluxes than the cuspier NFW, with the latter, as will be seen in what follows, 
being difficult to detect. The TN will be even more difficult to detect, implying for its 
detectability neutralino annihilation cross sections higher than the relevant range of values 
as predicted by supersymmetric models. 

Our final three choices represent well a range of possible concentrations of clumps in the 
galactic halo. In reality, the SIS configuration is very hard to achieve dynamically and has 



been recently ruled out for most of the clumps in our halo by EGRET bounds [ lOfl - We 
consider the SIS case to illustrate the range of clumps more highly concentrated than Moore 
et al. In addition, individual clumps may vary in concentration, given a range of possible 
initial conditions for a dark matter overdensity. Thus, an individual clump may be more 
highly concentrated than the mean. A possible range in concentrations should be kept in 
mind in general searches such as the one we study here. 

For a specific clump mass and distance, we need to determine the parameters of each of 
the three profiles. In the case of the SIS, we use the following two conditions to calculate p Q 
and the radius of the clump, R c i ump : 

1. The volume integral of the density must yield the mass of the clump: 

-R 

clump 

p(r) dV = M dump . (6) 
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2. The density of the clump at distance R c iump from its center must equal the local 
density, p G , of the galactic (host) halo at distance r c i ump , with r c i ump the distance of 
the clump from the center of the host halo: 

Pg(J clump) P clumpi^R clump) • (7) 

For the case of either the Moore et al. or the NFW density profile, we use Eqs.(|6|) and 
(0), but we need a third equation since we have three unknown quantities p Q , r s , and R c i ump . 
The third condition comes from the requirement that the radius of the clump be (smaller 
or) equal to the tidal radius (see, e.g., [32]). This guarantees that the clump is not tidally 



(8) 



stripped. Thus, 

r? _ p ^ r M c i urnp 

J^clump J^tidal — ' clump „ , / \ 

l3M G (r dump ) _ 

with M G {r dump ) the galactic halo mass inside a sphere of radius r dump . The SIS configuration 
is also stable against tidal stripping, as can be seen by considering the similar meaning of 
Eqs.flT]) and @. 

In both Eqs.(0) and (H) we need to specify the density profile of the galactic halo. Both the 
Moore et al. and the NFW profile are acceptable profiles for its description. Furthermore, 
for our purposes, they are essentially equivalent descriptions. Convergence studies have 
shown that these two density profiles are very similar at radii above ~ 1% of the virial 



radius (e.g., [0]). For a virial radius for our halo ~ 300 kpc, this means that the two 



density profiles are different only within the inner 3 kpc. Given that we are at 8.5 kpc from 
the galactic center, and given the nearest clump distances presented in Table |, the relevant 
distances from the center of the halo we are concerned with are far larger than 3kpc; thus, 
it is not of crucial importance whether the NFW or the Moore et al. density profile is used 
for the galactic halo . We choose the NFW profile to model the galactic halo, regardless 
of the density profile assumed for the clumps. For r s we use 27.7 kpc to match simulation 
results [f|. To find p c, the characteristic density for our galaxy, we normalize the density 
profile so that the peak circular velocity v is about 220 km/sec, and thus 



PoG 



ln(l + x Q ) - 



l + x 



(9) 



with x Q = R Q /r s , and R Q = 8.5 kpc our galactocentric distance. Finally, 



6.18 x lO 8 M /kpc 3 

Pair) — n- (10) 

V ; (r/27.7kpc) [1 + (r/27.7kpc)] 2 V ; 

One final characteristic of the dark matter clump structure is the core radius. Inside 
some radius R cor e, the annihilation rate becomes so large that the overdensity is destroyed 
as fast as it can fill the region, yielding thus a constant density core. Following ||, we find 
Rcore by equating the annihilation time scale to the free-fall time scale, 



T~ann Tff 



Py^core) i \ 

(av) 



2R 3 

" ' rl It 'Hip (11) 



GM, 



clump 



In this way we derive a maximum possible R core . The actual value might be smaller. Note 
that the choice of R core affects only the results for the SIS density profile since, as can be easily 
verified for the simple case of an unresolved clump, the flux coming from such a clump scales 
as Rcore ( see section [TV|) . In the case of the Moore et al. density profile the i? core -dependence 
of the flux is via the factor 0.67(l/x^ mp — 1) — 0.005 log x c i ump — logx core , approximately, 
with x c i ump = R c ium P / r s and x core = R cor e/ r s- Namely, the Moore et al. conclusions are not 
sensitive to the exact choice of R cor e |9 . In the case of the core-dominated fluxes of an SIS 
clump, using the maximum R COT e will yield the minimum fluxes. As will become clear from 
our results however, the conclusions do not change, at least qualitatively, since the SIS is 
steep enough to be detectable even when using these minimum fluxes. 

Below we discuss the annihilation cross sections into either continuum or monochromatic 
7-rays. We find that the dominant cross section times relative velocity product corresponds 



S 



to annihilation into the continuum, (<Jv) con t., and even though, in principle, all annihilation 
channels contribute in smearing out the central cusp, we use this dominant cross section in 
Eq.(|ITD to find R core . 

III. NEUTRALINO ANNIHILATION AS A GAMMA-RAY SOURCE 
A. The Neutralino-Supersymmetric models for dark matter 



Following |2y|, |35], |36|, [37| > we work in the frame of the Minimal Supersymmetric extension 
of the Standard Model (MSSM) using the computer code DarkSUSY p8| . The general 
R-parity conserving version of this model is characterized by more than a hundred free 
parameters. It is of common praxis to make some simplifying assumptions which leave only 
7 parameters, the higgsino mass parameter /i, the gaugino mass parameter M 2 , the ratio of 
the Higgs vacuum expectation values given by tan /3, the mass of the CP-odd Higgs boson 
ttla, the scalar mass parameter m , and the trilinear soft SUSY-breaking parameters A t 
and Ab for third generation quarks (for a discussion of the model and the simplification 



procedure, the results of which DarkSUSY incorporates, see |39], [40 



The LSP in most models is the lightest of the neutralinos. Neutralinos are a linear super- 
position of four neutral spin-1/2 Majorana particles: the superpartners of the Higgs, namely 
the neutral, CP-even higgsinos H\ and H2 , and the superpartners of the electroweak gauge 
bosons, namely the bino B and the wino W 3 . After electroweak symmetry breaking, these 
gauge eigenstates mix. The diagonalization of the corresponding tree-level mass matrix gives 
the neutralinos, 

= N ix B + N l2 W 3 + iV j3 ^i° + N l4 H 2 °. 

The lightest of these, the \ — Xi > is assumed to be the LSP and is referred to as the 
neutralino. Due to R-parity conservation the neutralino is stable, since there is no allowed 
state for it to decay into. Its mass is somewhere between some tens of GeV up to several 
TeV. It is cold, namely non-relativistic when decoupled, and is considered one of the most 
plausible candidates for the CDM particle. 

Using DarkSUSY, we made random scans of the parameter space, with overall limits for 
the seven MSSM parameters as given in Table 0. For the neutralino and chargino masses, 
one-loop corrections were used in accordance to |[42||, whereas for the Higgs bosons the 
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TABLE II: The ranges of parameter values used in the scans of the SUSY parameter space. 



Parameter 


H M 2 


tan/3 


m A m 


A b /m A t /m 


unit 


GeV GeV 


1 


GeV GeV 


1 1 


Min 


-50000 -50000 


3.0 


100 


-3 -3 


Max 


+50000 +50000 


60.0 


10000 30000 


+3 +3 



leading log two-loop radiative corrections were calculated using the code FeynHiggsFast |43 
incorporated in DarkSUSY. 

We scanned the parameter space to obtain the boundaries of the allowed SUSY parameter 
space which we plot in FIGs.1-3, 5, and 7. For each model, we checked whether it is excluded 
by accelerator constraints. We used the latest constraints available in the DarkSUSY pack- 
age, namely the 2000 constraints which are adequate for our purposes. The most important 
constraints come from the LEP collider at CERN, with respect to the lightest chargino and 
the lightest Higgs boson mass, as well as constraints from b — > 57 (for more details see, 



e.g., p5[ and references therein). For each model consistent with the accelerator constraints, 



we calculated the corresponding Q x h 2 , with Q x the neutralino relic density in units of the 
critical density, and h the Hubble parameter today in units of 100 km s -1 Mpc . Dark- 
SUSY calculates the relic density solving the Boltzmann equation numerically, while taking 
into account resonances and thresholds, and all the tree-level 2-body annihilation channels 

(see, 



p7| , j44||). We take into account the coannihilations between the lightest neutralino and 
the heavier neutralinos/charginos only if the mass difference is less than 30%, obtaining thus 
a relic density with accuracy ~ 5%, which is reasonable for our purposes. Lastly we are 
interested in models in which the corresponding relic abundance contributes a significant, 
but not excessive amount to the overall density. We choose the range 

0.1 < Q x h 2 < 0.3 (12) 

to be consistent with cosmological constraints. The lower limit in inequality (0) is chosen on 
the basis that for values lower than 0.1 there is not enough dark matter to play a significant 
role in structure formation, or to constitute a large fraction of the critical density. The upper 
bound is chosen on the basis that, given that there is a minimum age for the Universe, there 
is a maximum value that the dark matter density can have. The upper limit 0.3 corresponds 
to a minimum age of 12 Gyr approximately (see, e.g., pa]). Note that the choice of 0.3 
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for the upper limit can be considered generous, especially given the results from recent SN 
observations ^B, i7| , which constrain the allowed relic density to about Q x h 2 ~ 0.15, with 
h 2 ~ 1/2. 



B. Gamma-rays from neutralino annihilation 

Neutralinos annihilate through a variety of channels. We will focus on the continuum 
7-rays, and on the two monochromatic lines that are produced through the cascade decays 
of other primary annihilation products. The continuum contribution is mainly due to the 
decay of tt° mesons produced in jets from neutralino annihilations. Schematically, 

X + X — 7T° 

^7 + 7- (13) 

To model the fragmentation process and extract information on the number and energy 
spectrum of the 7-rays produced, we adopt a simplified version of the Hill spectrum [BR 



for the total hadron spectrum produced by quark fragmentation based on the leading-log 
approximation (LLA), 

\-"\i-x h y (14) 



dxh 16 

where Xh = Eh/m x , and is the energy of a hadron in a jet with total energy equal to 
the neutralino mass, m x . Assuming that all the hadrons produced are pions, which is very 
close to reality, and assuming that each pion family takes approximately 1/3 of the total 
pion content of each jet we may write for the 7r° spectrum, 

i£ (15) 

with x n = E n o/m x . Furthermore, the probability per unit energy that a neutral pion with 
energy E^o produces a photon with energy E 1 through the process shown in Eq.(|TB|) is 
2/E n o. Thus, from Eq.fllSD we get for the continuum photon spectrum, 

dN cont . f 1 2 dN„o 



/ --r-dy 16 

Jx-, y dy 



with x 7 = E 1 /m x and y = E v o/ 



m 



The continuum signal lacks distinctive features; this might make it difficult to discriminate 
from other possible 7-ray sources. A more unique signature is given by monochromatic 7- 
rays, which arise from loop-induced s-wave annihilations into the 77 and the Zj final states. 
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TABLE III: VERITAS specifications 



Energy range 50 GeV-50 TeV 

Energy resolution 15% 

Effective area Angular resolution 

1 x 10 8 cm 2 5' at 100 GeV 

4 x 10 8 cm 2 3' at 300 GeV 

1 x 10 9 cm 2 2' at 1 TeV 



These lines are free of astrophysical backgrounds. The respective spectra for the two lines 
are: 

• For the 77, 2 photons at ~ m x 

• For the Z7, 1 photon at E 7 ~ m x — m 2 z /4:m x , 

as can be easily verified by the conservation laws and the fact that the neutralinos are 
expected to be highly non-relativistic. 

Note that for high m x values, the energies of the two lines are expected to coincide. 
Furthermore, in the case of the Zj line, the energy of the photon becomes vanishingly small 
as m x tends to rriz/2. Below m x = mz/2, the Z7 process is no longer kinematically allowed. 



IV. THE GAMMA-RAY SIGNAL 



A. Fluxes and counts 



We derive the fluxes for resolved and unresolved clumps as a function of the neutralino 
parameters, and the size and distance of the clumps, assuming the nearest clump for each 
mass bin. To turn the estimated fluxes into photon counts collected by an ACT, we need 
to specify the characteristics of the particular ACT. We choose to use VERITAS as the 



standard next generation ACT and list its specifications in Table |TJ 
The number, R, of annihilations per unit volume and time is 



R = n 2( av ). (17) 
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where n is the number density of neutralinos, and (av)i is the thermally averaged product 
of the cross section times the relative velocity ]5(J for the continuum 7-rays (i=cont.), the 
27 line (2=77), and the Z7 line (i = Zj). Denoting by iVj the number of 7-rays emitted per 
annihilation, with energy E 1 above the energy threshold E t h used, then iV 77 = 2 (provided 
that m x > E t h), and Nz-y = 1 (provided that m x — m 2 z /4m x > Eth)', N cont . can be obtained 
by integrating the spectrum given by Eq.fllED, 

N cont _{E 1 > E th ) = [ dNcont - rf x = ^ x 3/2 _ ^ + 5^^- + ^ — - ^ . (18) 

Jxth dx i 6 3 3 

where x th = E th /m x . 

The emission coefficient ji, namely the number of emitted photons per unit time and 
volume, will be 

A = N<*{ov) t = (19) 

where p is the mass density of neutralinos. 

Given the expression for the emission coefficient, we can calculate the flux coming from 
a clump. The observed flux at the Earth depends on whether the object is resolved or 
unresolved by the ACT, namely, whether the angular size of the clump is larger or smaller 
than the angular resolution of the instrument. More specifically, in the case where the clump 
under consideration is unresolved, its flux at the Earth, F^ nres , is given by the volume integral 
of the emission coefficient divided by the area Arcd 2 , 

= 55*/ ^* = ^nsjr/ ^ R (20 » 

where d denotes the distance of the clump from the Earth and R c iump is, as before, the clump 
radius. At this point, note that d and r dump (appearing in Eqs.(|7|) and (§)) are two different 
quantities; d is the distance of the clump from the Earth, whereas r c i ump is the distance of 
the clump from the center of the Milky Way. For a clump with galactic coordinates (/,&), 
the two quantities are related via the expression, 

r dum P = d 2 + R 2 Q - 2dR Q cos I cos b (21) 

with R =8.5 kpc. If the clumps are resolved, the corresponding flux at the Earth, F* es , is 
given by the integral along the line of sight of the emission coefficient, 

FL = i / j,dz = f <Uz (22) 



1-i ./;„. ' lr; in 2 
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In the case of the unresolved sources, the source counts, Nl nres , will be simply, 

^unres = F unres A e fft (23) 

with A e ff the effective collecting area of the instrument and t the integration time. For a 
resolved clump, we can define its surface brightness /i* as seen by an ACT as 

A** = KesAefft (24) 

where A e ff and t the same as above. To obtain the number of photons, N* e8 , that the 
instrument will collect from the source, assuming that we observe the central pixel, we 
integrate the surface brightness of the source over a disk centered on the source with angular 
radius equal to the angular resolution of the instrument <j e , 

i-2-k PCTg 

K es = / d4> / ftOM. (25) 



Jo JO 

The angular resolution of the instrument defines a maximum projected size R ma x for a clump 



at distance d by ag = R max /d. Similarly, 6 = R/d. Substituting Eqs.(|22|) and fl2"4" ) in Eq.(p5| 
we obtain, 

1 Niiav) 



K es = — 2 ^A eff tJ G d^j Q RdRjJdz. (26) 
We find that most clumps we considered will be resolved by future ACTs, e.g., VERI- 



TAS f5"Ifl . Thus, we focus on the resolved fluxes and counts when discussing our results in 



^V. For the diffuse 7-ray contribution of clumps in the halo of our galaxy see & [H| 



B. Background counts 

The detectability of the clumps depends on their fluxes as compared to the possible 
background contributions. As we discuss below, the dominant background contributions are 
due to electronic and hadronic cosmic ray showers. Efforts to discriminate between photon 
primaries versus charged primaries can help improve considerably the detectability of dark 
matter clumps. We also discuss the galactic and extragalactic 7-ray backgrounds, as well as 
the contribution from the rest of the dark matter halo component. 
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1. Electronic and hadronic cosmic ray showers. 



For the background rate of 7- like hadronic showers we use the integral spectrum presented 

in 



52 



-1.7 



^>«^.l,r=(A) -W. (27) 

This was derived from data taken with the Whipple 10m telescope. VERITAS is expected to 
have an improved method for rejecting charged primary showers, compared to the Whipple 
telescope. However, aiming at a conservative treatment, we do not take this into account. 



For the electron-induced background we use the following integrated spectrum |53 

The electron-induced showers at Whipple are indistinguishable from 7-rays and can be 
rejected only on the basis of their arrival direction, whereas the hadron-induced showers 
are more extended on the ground than the electron-induced ones. Furthermore, we see 
that cosmic-ray electrons have a steeper spectrum than cosmic-ray nuclei and dominate the 
background at low energies. 



2. Extragalactic and galactic ■j-ray emission. 

For the extragalactic diffuse emission, we use a fit to the EGRET data, valid for the 
energy range from 30 MeV to ~ 100 GeV p>3 



rIN / F \ -2.10±0.03 

= (7.32 ± 0.34) x 10- 9 — %- - MeV-W W- 1 . (29) 

dVLdE y ' V 451MeV / 

The galactic diffuse emission is thought to be mainly due to cosmic-ray protons and 
electrons interacting with the interstellar medium. It is enhanced towards the galactic 
center and the galactic disk, and has been measured by EGRET |53], [K| up to about 10 
GeV. For its differential spectrum we use the expression provided in J52"| |, namely a power 
law fall-off in energy of the form 

^r(E,l,b) = 10" 6 iVo(/,&) GeV-W-'s" 1 *- 1 . (30) 

The factor N depends only on the galactic coordinates (/, b) and has been fixed using 
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EGRET data at 1 GeV (see [52]], and references therein). It is given by 



N (l, b) = - 85,5 + 0.5 if\l\> 30° 

v/l + (1/35) Vl + (6/(1-1 + KI0.022)) 2 " 

+ 0.5 z/ |/| < 30° (31) 



v/l + (//35)Vl + (6/l-8) 2 
with / and b varying in [—180°, 180°] and [—90°, 90°], respectively. This expression is in 
reasonable agreement with data especially towards the galactic center. The exponent a in 



principle depends on both the energy range and the galactic coordinates; as in |52| , p6[ , we 
use a = —2.7, assuming the behavior continues at higher energies. 

It is also important to obtain an estimate of the 7-ray counts due to neutralino annihi- 
lation in the smooth dark matter component of the galactic halo. The flux of continuum 
7-rays of the smooth halo component is given by 



1 N cont ( crv) 
47r mt 



^smoothW = — C ° m 2 N ' I PhmW- (32) 



with pc given by Eq. (|T0D , and if) the angle of the direction of observation with respect to 
the direction of the galactic center. This flux is clearly enhanced towards the central regions 
of the galaxy (for more details see, e.g., 0). In order to estimate a value close to the 
maximum possible smooth component counts, N smoot h, we use the direction ip = 10° - close 
to the center. The flux coming from this direction is approximately 



F S mooth(E > E th ) ~ 3.5 x \0 N con i ,con ^ s^cnT^r" 1 . (33) 

The emission due to the total clumped component of the dark matter halo, as seen in 



|9|, overwhelms the smooth component 7-ray emission. Note that neither the smooth, 
nor the clumped component 7-ray emission constitutes an additional background for our 
study. For example, part or all of the extragalactic background, which we have already 
taken into account, would be due to numerous, unresolved clumps in our halo. In fact, this 
has allowed ref. [10| to place strong limits on the neutralino parameter space for clumped 



halos with SIS and Moore et al. clump profiles, by comparing the clumped component 7-ray 
emission with the extragalactic diffuse 7-ray emission (Eq. (|29|) ). In particular, an SIS clump 
distribution would surpass EGRET bounds for all of the SUSY models usually considered 
for neutralinos. 

In order to estimate the relative importance of the several contributions we assume an 
ACT observation with the following characteristics: ^4 e // = 10 8 cm 2 , t = 100 hrs, E t h = 50 
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GeV, ae = 0.1°. For these parameters, we find the counts iV e _ ~ 1278, N h ~ 2718, N eg ~ 6, 
Nq — 225 (for m x = 100 GeV and, referring to the continuum 7-rays, for a typical cross 
section (<rv) cont . = 5 x 10 _27 cm 3 s _1 we obtain N smooth ~ 0, rounding off to the closest 
integer). To calculate the galactic background counts, we made the conservative assumption 
of the maximum possible background and, thus, used (7,6) = (0°,0°). This background is 
anisotropic and will be lower at higher galactic latitudes. Evidently, the main backgrounds 
for ACT searches are the hadronic and electronic backgrounds. Only if the composition 
of primaries can be identified, the galactic and the extragalactic diffuse 7-ray backgrounds 
become important; this would be highly desirable, given that the galactic background is (at 
least) an order of magnitude lower than the cosmic ray induced backgrounds. 

We chose to keep a conservative estimate of the backgrounds by considering the two 
dominant backgrounds, the electronic and the hadronic cosmic-ray induced ones, since, as 
we have shown above, the other contributions are at least one order of magnitude lower. 
Thus, the background counts in the case of the continuum are 

N b (E > E th ) = ^(E > E th )AQ A eff t (34) 

where we define dNb/dfl as the sum of the hadronic and the electronic cosmic ray induced 
backgrounds, as given by Eq.(^J) and (|28"D, respectively. In the case of the lines the back- 
ground counts are given by 



N b (E = E ) 



dNh /ri „ . dNi, . „ . 
(E > Ei) — -jtt(E > E 2 ) 



Att Aeff t (35) 



dn y ~ ' dtt 

where E a stands for the energy of the line and Ei, E 2 stand for (1 — R E )E a and (1 + Re)E q 
respectively, with R E denoting the energy resolution of the detector. 

Note that the only anisotropic background - the galactic - is negligible in comparison 
with the charged particle backgrounds. This makes our calculations essentially independent 
of the exact direction of the clump in the sky. Any particular choice of galactic coordinates 
would not limit the generality of our conclusions with respect to the observability of the 
clumps. The only assumption we should make is that the clumps are more likely to be 
found away from the galactic plane. This is to be expected since simulations show that the 
orbits of CDM satellites in present-day halos very rarely take them near the disk (e.g., [f22||). 
This is not a very constraining assumption - all directions in the sky that avoid the galactic 
disk are essentially equivalent - and thus, the only really important quantity with respect 
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to the location of a clump is its distance from us. The exact galactic coordinates can in 
principle be important via Eqs.(^) and fl2~T|), but it turns out that this is not the case. For 
example, considering the clumps with mass M dump = 10 8 M Q and nearest distance of about 
4 kpc (see Table |), we find that putting such a clump towards the galactic center yields 
p G ~ 2.8 x 10 9 M Q kpc -3 , whereas putting it in the exact opposite direction, towards the 
anticenter, yields pa — 1.3 x 10 9 M Q kpc -3 , namely there is no significant variation. 



V. RESULTS 



A. The detectability condition 

The background follows Poisson statistics and as such it exhibits fluctuations of amplitude 
y/Nb- For given technical and observation parameters, the minimum detectable flux of 7- 
rays for an ACT is determined by the condition that the significance S exceeds a specific 
number, M s , of standard deviations, a, or that the number of detected photons exceeds a 
specific number N Q . Namely, 

N 

S = > M s (36) 



and, 

N s > N . (37) 

The first condition simply defines what achieving an M s -a detection requires. The second 
condition ensures that we obtain sufficient statistics of 7-ray events and is useful especially 
at the higher energies ( > 1 TeV), where the backgrounds are negligible, and thus the first 
condition is useless. The second condition also shows that at high energies, whether the 
source will be detectable or not is determined only from its high energy spectrum (and the 
ACT characteristics). Given that we start at tens of GeV energy thresholds, we check for 
the first condition making sure that the second condition is satisfied as well. We require a 
5-0" detection level (M s = 5). As has been already mentioned, for all the mass bins studied, 
the nearest clump can be resolved for standard ACT angular resolutions (see Table [IT]] ), and 
thus we use Eq.(^) for N s , whereas for the background we use either Eq. (|3~4D or Eq.([35|), 
depending on whether we study the continuum or the lines, respectively. 

Using the aforementioned detectability condition, we can constrain the SUSY parameter 
space. In addition to the unknown neutralino parameters, (crv)i and m x , we do not know 
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the precise density profiles of the dark matter clumps. In the figures that will follow, we take 
under consideration the profile uncertainties by using the three density profiles mentioned 
before (the SIS, Moore et al., and NFW profiles). We then set the constraints on the 
neutralino parameter space as follows: we substitute Eq.(^) for N s , and Eq.(Q) or ( |3"5] ) 



for Nf,, as well as the expressions for iV, in condition (^), and then we solve the resulting 
inequality with respect to {av)i. This yields an inequality of the form 

(<ro)i > /(/, A eff , t, Efh, dN b /dn, An, M s )m 2 x . (38) 

In other words, an M s -a detection for given A e ff, t, E t h, dN^/dQ, AQ, and M s can be 
obtained as long as (crv)i is above a certain value which depends on m 2 and the integral of 
the density squared, which we have denoted by I. Thus, plotting inequality (|38| ) with the 
equality sign onto the SUSY parameter space, we divide the SUSY parameter space into the 
detectable (above the line defined by (|38[) ) and the undetectable (below the line defined by 
(PS])) regions. 

In another set of figures, we show the flux for a given clump as a function of m x , assuming 
a typical value, (av) typ , for the cross-section. This is interesting to be compared to the 
dominating backgrounds, as well as to what we call the minimum detectable flux F m i n , 
which we define using again the detectability criterion given in (^). More specifically, for 
the minimum required source counts N s ^ min , so that an M s -a detection level be achieved we 
will have from fl36|) 

N s , min = M s ^/N~ h (39) 

For N s>min = F min A eff t we obtain, 

F • = (40) 

± mm t ^*±vj j 

A efft 

This is separately applicable for the continuum and the lines. 



B. Continuum 7-rays 

The minimum detectable neutralino annihilation cross section into continuum 7-rays as a 
function of neutralino mass for a 10 2 M Q and a 10 8 M Q dark matter clump is shown in FIG.[l|. 
These curves were obtained using condition (|38|) for A e ff = 10 8 cm 2 , t = 3.5 x 10 5 s (100 hrs), 
E th = 50 GeV and ag = 6' (these are also the parameters we use for the rest of the figures, 
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unless otherwise stated). In the same figure, the three cases of clump modeling using an 
NFW, a Moore et al., or an SIS density profile are shown. Note that by plotting the results 
for the minimum and the maximum clump masses that we considered, 10 2 M Q and 1O 8 M , 
respectively, any other clump mass choice lies in between. The same is true for the density 
profiles, at least with respect to the degree of central concentration: any profile steeper than 
the NFW and softer than the SIS will be a curve between the two limiting curves defined 
for a specific clump mass, by the SIS and the NFW profiles, and as an example we plot the 
Moore et al. case. 

In the case of a 10 2 M Q clump, we see that if the clump were described by a density 
profile as centrally concentrated as the SIS, then it would be detectable for most neutralino 
parameters, unless m x '^ J 70 GeV. The same clump modeled with the Moore et al. profile has 
good chances to be detected, since for a large range of m x the (<Jv) con t. — rn x curve is below 



most of the possible models. Furthermore, note that given the recent results from [10], these 
highly centrally concentrated clumps would have to be a rare event among concentrations 
of halo clumps. 

The same results are plotted for a clump with mass 1O 8 M . The minimum detectable 
(<7v) cont . for all the three profiles goes down by approximately two orders of magnitude, 
compared to the 1O 2 M case. As a result, such a massive clump is most likely detectable, 
even if it is best described by the least concentrated density profile, the NFW. Note also 
that, referring to the low mass clumps, the NFW case is, compared to the other density 
profiles, the most consistent profile with EGRET bounds . 

The way the constraints imposed on the SUSY parameter space depend on the mass of 
the clump is depicted in FIG.^], assuming a Moore et al. density profile, and for the same 
observation parameters as in FIG.[l]. The results are plotted for four different clump masses: 
10 2 M Q , 1O 4 M , 10 6 M q and 10 8 M Q . As evident, the higher mass clumps have better chances 
of detectability, since the higher the mass of the clump, the larger the part of the SUSY 
parameter space that lies above the minimum detectable cross section, ((Jv) cont , versus m x 
curve; equivalently, the higher the mass, the stronger the constraints on the SUSY parameter 
space that the non-detectability of an otherwise detected clump (via e.g. its synchrotron 
emission |Tl|]) imposes, given that if the clump is not detected by ACTs, all the models 
above the {crv) cont .- r m x curve will be excluded. 

The minimum neutralino mass that can be explored also depends on the observation 
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features, and more importantly, on the energy threshold. To show the role of the energy 
threshold, E t h, in the detectability of the clumps for the SUSY parameter space, in FIG.|3| 
we plot the same kind of curves as in previous figures, using a 10 5 M Q Moore et al. clump 
and for three different energy thresholds: 50, 100 and 250 GeV, for 100 hrs of observation 
and for A e ff and <jq in agreement with Table [TIT]. Using a higher energy threshold means 
both a larger A e ff and a better angular resolution, as well as lower contributions from the 
backgrounds, which makes it clear that ACTs are well-suited to explore the higher m x range 
of the neutralino parameter space. On the other hand, the lower the energy threshold the 
larger the part of the parameter space that ACTs can explore, but even in the best case of the 
smallest energy threshold, the low mass parameter space which corresponds to a broad range 
of (av) cont . will remain unexplored regardless of the mass of the clump. This demonstrates 
the need for complementary observations at lower energies, such as observations with GLAST 



57], accelerator searches, and direct searches. 



Our second set of results concerns the flux coming from a clump as a function of neutralino 
mass. Following the trend given by the DarkSUSY models, we choose a typical cross section 
for the annihilation into the continuum (<Jv)t yp — 5 x 10 -27 cm 3 s . Note that the SUSY 
parameter space is not extremely wide in the continuum case, at least for neutralino masses 
above ~ 80 GeV, and thus the choice of a typical cross section is reasonable. The flux as 
a function of m x for (av) typ , and for clump masses M dump = 1O 2 M and M c i ump = 1O 8 M 
is shown in FIG.|j. As above, we show the results for the three different density profiles we 
used. The energy threshold was set to 50 GeV. On the same plot and for the same energy 
threshold, the electronic and hadronic backgrounds, as well as the 5-cr minimum detectable 
flux, F min (defined in Eq.(|(J)), are plotted using A eff = 10 8 cm 2 and t = 3.5 x 10 5 s (100 
hrs). In these plots a clump of a specific mass and density profile will be detectable for those 
neutralino masses for which its flux is greater or equal to the minimum detectable flux F min . 
We will call this range of neutralino masses the m x detectability range. 

From FIG.[|, we see that a 1O 2 M NFW clump will not be detectable (regardless of 
m x ), an SIS clump of the same mass has an m x detectability range m x ^55 GeV, whereas 
if the clump density distribution is better described by a Moore et al. profile, the m x 
detectability range will be 87 GeV ^,m x ^2.9 TeV, approximately. In the case of a 1O 8 M 
clump, the SIS profile has an m x detectability range m x ^53 GeV. Similarly for the Moore 
et al. 1O 8 M clump m x ^65 GeV, whereas the same clump if modeled with an NFW profile 
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will have an m x detectability range 77 GeV ^m x ^5 TeV. Focusing on the intermediate 
case, the Moore et al. profile, we see that in terms of constraints on m x , the detectability 
of low mass clumps provide both a lower and an upper limit (i.e., detectability range of 
the form Tn X)min %m x %m Xjmax ) , whereas the higher-mass clump detectability provides only 
a weaker lower limit (i.e., detectability range of the form m x ^m Xtmin /). Note that the 
same conclusions if applied in the case of a clump known to exist but which turns out to 
be non-detectable by ACTs, imply complementary to the above limits for the neutralino 
parameter space. Namely, the non-detectability of the clump will imply that m x lies in 
the range : all m x values predicted by SUSY — detectability range. Actually, as can be 
seen from all figures referring to the continuum signal, it will be the non-detectability that 
will impose the strongest constraints on the SUSY parameter space. Furthermore, although 
the larger the clump mass the easier to detect, having a range of clump masses is needed 
in order to determine the neutralino parameters. As we discuss below, the constraints 
can become stronger when information on the line emission is combined with the better 
detectable continuum 7-rays. 

Lastly, given the SUSY parameter space and the (<Jv) cont ./m x dependence of the flux, we 
find the minimum F min and maximum F max flux above 50 GeV that a clump of a specific 
mass and density profile can emit. We do that by finding from the SUSY parameter space 
what {av) cont . and m x combinations yield the minimum and maximum (crv) con t./m x ratio, 
respectively. In this way, we bracket the flux of a clump and we know that, even for density 
profiles other than the three ones used here, as long as they are somewhere between the 
NFW and the SIS, their flux will be somewhere in the interval [F min , F max }. The results 



of this calculation are shown in Table [TV]. For the minimum possible clump flux, F min , we 
chose m x = 5 TeV and (av) = 10~ 27 cm 3 s _1 . This combination yields a representative 
lower flux value, at least for the horizontal feature of SUSY models in the large m x range 
(even lower fluxes are obtained for the low mass and low cross section part of the SUSY 
parameter space, but this range is not likely to be accessible by ACTs, anyway). For the 
maximum possible clump flux, F max , we used m x = 60 GeV and (av) = 2 x 10~ 26 cm 3 s _1 , 
which is a combination that quite reliably gives the maximum possible flux. For comparison, 
in Table [V] we give the minimum required flux, F miritdet , so that a 5-cr detection be obtained 
for certain E t h, A e ff and t combinations, as was calculated using Eq. (EH) and assuming the 
hadronic and electronic backgrounds as the dominant ones. 
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TABLE IV: Minimum and maximum possible fluxes for the three different density profiles - see 
text for details. All fluxes are in photons cm _2 s _1 and masses are in M & . 
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TABLE V: Minimum detectable fluxes for specific energy thresholds, effective areas and integration 
times so that a 5-cr detection be obtained. All fluxes are in photons cm -2 s -1 , energies in GeV, 
areas in cm 2 and integration times in s. 

Eth A e ff t Fmin,det. 

50 1.0 x 10 8 3.5 x 10 5 8.9 x lO" 12 

100 1.0 x 10 8 3.5 x 10 5 4.6 x lO" 12 

250 1.0 x 10 s 3.5 x 10 5 2.0 x 10- 12 

50 1.0 x 10 9 3.5 x 10 5 2.8 x 10- 12 

50 1.0 x 10 8 2.6 x 10 6 3.2 x lO" 12 

50 1.0 x 10 8 3.2 x 10 7 1.1 x 10- 13 



From the two tables we see that the least likely case of clumps, namely clumps as dense as 
the SIS, can be easily detectable in less than 100 hrs. In the case of the Moore et al. profile, 
the higher mass clumps (> 1O 6 M ) will be easily detectable since both F min and F max exceed 
the minimum required flux, F minjdet _, to achieve a 5-cr detection level. The lower mass clumps 
have an easily detectable F max , but not necessarily a detectable F min . A 10 4 M o Moore et 
al. clump would have a detectable F min for either a large integration times (2.6 x 10 6 s ~1 
month ), or for larger collective areas (e.g., 10 9 cm 2 at 50 GeV), which is not among the 
abilities of both existing and upcoming ACTs. Finally, NFW clumps with M c i ump ^10 6 M & 
are detectable and this is particularly important given the EGRET constraints placed on 
the most concentrated clumps, such as SIS and Moore et al. ones. In contrast, even a whole 
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year of observation (3.2 x 10 7 s) would not be enough to detect smaller mass NFW clumps 
with fluxes that are closer to the corresponding F min . 

C. Monochromatic 7-rays 

The minimum detectable cross section (cru) 77 for annihilation into monochromatic 7- 
rays with energy equal to m x , as a function of m x is shown for M c i ump = 1O 2 M and 
M c i ump = 1O 8 M in FIG.|5|. The constraining curves stop at E t h, which in this case is 
50 GeV (below this the cross section is non-zero, but the number of photons with energy 
above 50 GeV coming from neutralinos with mass less than 50 GeV is, of course, 0). A 
first, clear difference, in comparison with the continuum 7-rays, is the fact that the SUSY 
parameter space is significantly wider with respect to the cross section. Therefore, imposing 
constraints on the SUSY parameter space using the detectability (or non-detectability) of 
the clumps through their 77 line emission will be harder than through their continuum 
emission. Nevertheless, if lines are detected in addition to the continuum, the degeneracy 
in the SUSY parameter space will be lifted since a narrower range of models would fit both 
detections. 

This is depicted more quantitatively in FIG.|| where we present the flux of monochromatic 
7-rays as a function of m x , as well as the minimum detectable flux and the two major 
background contributions. Note that in this case both the minimum detectable flux and 
the background fluxes depend on m x , given that the energy of the monochromatic 7-rays is 
equal to m x . As we did in the case of the continuum, to find the flux coming from a given 
clump we assume a typical cross section, in order to express the flux only as a function of 
m x . Based on the SUSY parameter space presented in FIG.[|, we choose the large cross 
section region and use (crv) 77 = 2 x 10~ 30 cm 3 s _1 . This value yields a flux close to the 
maximum possible flux such that a comparison can be attempted with the continuum case, 
which is generally stronger. We see, however, from FIG.|| that even in this case, the 77 
line will not be observable with upcoming ACTs. More specifically, for M c i ump = 1O 2 M the 
chances of observing such a light clump in the monochromatic 7-ray lines are bad, regardless 
of neutralino mass. With respect to the density profile, neither the Moore et al. nor the 
NFW will be detectable. In the case of a density profile between the Moore et al. and the 
SIS, observability could be achieved imposing at the same time an upper limit on m x . In 
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the case of M c i ump = 1O 8 M , a Moore et al. clump will be detectable as long as m x ^ 4.5 
TeV, whereas the same clump, if described by an NFW profile, would be detectable only if 



m x <68 GeV. 



The case is similar for the Zj line. Results for the SUSY parameter space are presented 
in FIG.[7|, for M c i ump = 10 2 M Q and 1O 8 M . In these figures we require that > E th . Using 
the fact that E 1 = m x — m z /4m x , this implies that the neutralino mass 



E th + V m l + E , 



m x > ^ ^ (41) 

which translates to m x ^77 GeV for E t h = 50 GeV. The two lines are expected to coincide 
at relatively large m x , and thus, the conclusions for both are similar with the difference that 
due to the above constraint the Zj line explores the parameter space starting from higher 
neutralino masses compared to the 77 line. In addition, the SUSY models corresponding to 
m x above 500 GeV are less dispersed in the case of the Z7 line than in the case of the 77. 
Again, this will help lift degeneracies of neutralino models, if any clumps are detected. 

The flux as a function of neutralino mass assuming a maximum (crv)z-y = 3 x 10~ 30 cm 3 
s _1 , is shown in FIG|| In general the conclusions are the same as in the 77 line. The light 
clumps are not very promising with respect to their detectability. A 1O 8 M Moore et al. 
clump will be detectable as long as m x ^, 6 TeV. Note, however, that this is true under the 
assumption of the maximum (crv)z-y, and thus of the maximum possible flux. 



VI. DISCUSSION AND CONCLUSIONS 



The results presented above were based on specific choices for the detector and the ob- 
servation parameters. It is worth noting that these results can be easily scaled with some of 
these parameters. For instance, the effective area A e ff and observation time t have a simple 
scaling. The minimum detectable (av) for a specific m x scales as l/(A e fft) 1 ^ 2 . An order 
of magnitude lower (erf) could be explored using two orders of magnitude higher A e fft, for 
example, A e ff = 10 9 cm 2 and t = 3.5 x 10 6 s (~ 40 days). This would be very helpful in 
the case of the hard to detect lines (see FIG.| and FIG.0). If some clumps are detected in 
the continuum, a reasonable strategy would be to subsequently search for the weaker line 
signal over a longer t. Similarly, the minimum detectable flux plotted in FIGs.^, and |8|, 
depends on the combination A e fft. With respect to the angular resolution erg, the minimum 
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detectable (av) has a dependance of the form oc l/yAfi, approximately, while the mini- 
mum detectable fluxes are oc V Af2. Note that although the nearest clumps are resolved, 
we used the angular resolution of the instrument and not the angular size of the clumps 
as the relevant angle for the collection of background noise, assuming that we observe the 
central pixel of the source - where most of the source counts originate from - treating thus, 
the source with respect to the backgrounds as being effectively unresolved. If, instead, one 
uses the angular size of the object as the relevant angle for the collection of noise, both the 
detectable cross sections and the minimum detectable fluxes are about two and one orders 
of magnitude higher for a 10 8 M Q and a 10 2 M Q clump, respectively. 

An essential parameter in our calculations is the distance of the source from the Earth, 
which we assumed to be the estimated distance to the nearest clump. However, our results 
can be easily extended to larger distances or, equivalently, to lower numbers of clumps given 
the l/d 2 dependence of the flux, which is actually the most significant dependence of our 
results on the distance [58|]. For example, one order of magnitude larger distance makes the 
minimum (av) that can be explored two orders of magnitude higher; the same is true for the 
minimum detectable fluxes. This is very useful given the recent bounds on SIS and Moore 



et al. clumps from EGRET data ||10|| . More specifically, only NFW clumps can populate 
the whole halo as seen in CDM simulations without overwhelming EGRET limits. For 
the other density profiles, consistency with EGRET data could imply, e.g., a smaller than 
currently believed clump abundance, which would translate, in our case, to larger distances 
of the clumps. Using the above scaling of our results with distance, it is easy to conclude 
with respect to the clump detectability in any other case that would change the distance 
estimate. Similar is the case for the density profile assumed for the clumps; the detectability 
of a density profile, e.g., with central concentration in between that of the Moore at al. and 
the NFW profiles, can be predicted, at least qualitatively, using our results. 

The choice of the energy threshold will be a crucial factor determining what will be the 
minimum m x values accessible to ACTs. An energy threshold low enough, along with the 
ability of ACTs to identify primaries leading thus to backgrounds smaller by an order of 
magnitude, are highly desirable for future ACTs, if they are to be used to explore the SUSY 
parameter space and detect even the less concentrated clumps such as the NFW case. 

In addition, we would like to re-emphasize the strategy of combining observations in the 
continuum and the annihilation lines for different clump masses in order to constrain more 
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effectively the SUSY parameter space. For example, from FIGs|4| and |6| we saw that a low 
mass clump has an m x detectability range of the form m m i n % i m x % J m max . Combining this 
with line observations, one could narrow the range of neutralino masses. More precisely, one 
should study the allowed models in SUSY parameter space that fit the detected continuum 
flux and check which ones give a detectable line signal. 

The detectability of the lines can yield m x detectability ranges of the form m x < m' max 
for various combinations of clumps masses and density profiles, assuming a typical or an 
optimum (erf) for the neutralino annihilation. Assuming for the lines the maximum (crv), a 
10 8 M© Moore et al. clump will be detectable both in the continuum and in the lines, if m x 
is somewhere between 65 GeV and 4.5 TeV, where rn x ^ 65 GeV comes from the continuum 
limit and m x ^4.5 TeV comes from the lines. Although not a very narrow interval in m x , 
excluding the very low mass end of the allowed neutralino mass range in combination with 
an upper m x -limit for detectability, can be very useful. It would narrow down the range in 
(ffv) cont . significantly. Another example is that of a 10 8 M Q NFW clump which is detectable 
in either the continuum ( detectability range: 77 GeV < m x ^ 5 TeV) or the 77 line 
(detectability range: m x < 68 GeV), but not both. This fact can, in principle, be used to 
extract information about the mass of the clump. 

Our study of the detectability of clumps by ACTs is closely tied to the lack of knowledge 
with respect to the SUSY parameters. Thus our results can be used in two ways, depending 
mainly on what will be determined first, the neutralino parameters or the distribution and 
profiles of CDM clumps. If clumps are discovered first (for example, through synchrotron 
emission or from GLAST 7-ray observations), 7-ray studies by ACTs can help narrow the 
neutralino parameter space. In most cases, the strongest constraints on the SUSY parameter 



space arise from the non-detectability in 7-rays of a clump seen, e.g., in the microwaves fll 
rather than its detectability as, for example, in the case of a 10 8 M Q Moore et al. clump. 
Such a clump will be detectable in both the continuum and the lines if m x is between 65 and 
4.5 TeV. This large range of neutralino masses can be ruled out if the clump is not detected 
by VERITAS. This can become quite useful in the future. For the time being, what is even 
more important is that such a clump, if existent, will be easily detectable almost for any 
possible neutralino mass. 

Summarizing, we conclude that the chances of detectability of dark matter clumps in the 
galactic halo, due to continuum 7-rays produced by neutralino annihilation via upcoming 
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ACTs are very good, especially in the case of relatively massive and highly centrally con- 
centrated clumps. The signatures expected in the form of the 77 and the Z7 lines are less 
easily detectable, but can be quite useful in lifting degeneracies among neutralino SUSY 
parameters. If we figure out the neutralino parameters first, ACT searches will be strong 
tests of the CDM model predictions with respect both the substructure in galactic halos and 
the density profiles of the substructure clumps. ACT searches will be able to explore large 
ranges of neutralino masses above ^ 50 GeV (or more, depending on the energy threshold). 
These studies will complement the searches for lower mass neutralinos, accessible through 
high energy accelerators, direct and other indirect detection methods, and satellite 7-ray 
telescopes such as GLAST. 
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FIG. 1: The minimum detectable (crv} con t. versus m x for the SIS, the Moore et al., and the NFW 
profile. The clump masses used are lCP-M© and 1O 8 M0. The dots represent allowed SUSY models 
(see text for details). The lines represent the 5-cr detection for A e ff = 10 8 cm 2 , E t h=50 GeV, 
ag = Q' , and for 100 hrs of observation. Only SUSY models that lie above the corresponding curve 
will yield a detectable signal. 
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FIG. 4: The flux of continuum 7-rays as a function of m x , for (crv} con t. = 5 x 10 27 cm 3 s 1 , and 
for the SIS, Moore et al., and NFW profile. The clump masses used are 1O 2 M and 1O 8 M . Also 
shown are the two dominant background contributions, the hadronic cosmic ray and the electronic 
cosmic ray induced ones. The dashed line represents the minimum flux of the clump at the Earth, 
required so that a 5-cr detection level be achieved for A e ff = 10 8 cm 2 , E t h=50 GeV, a$ = 6', and 
100 hrs of observation. For m x values corresponding to fluxes higher than F m i n , the clumps will 
be detectable in the continuum. 
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FIG. 5: The minimum detectable (av)^ as a function of m x for the SIS, the Moore et al., and 
the NFW profile. The clump masses used are 1O 2 M0 and 1O 8 M0. The dots represent allowed 
SUSY models. We require a 5-a detection level for A e ff = 10 8 cm 2 , £^=50 GeV, oq = 6 , energy 
resolution Re = 15%, and for 100 hrs of observation. For a specific clump mass and a specific 
density profile, and for the above instrument and observation characteristics, only the SUSY models 
that lie above the corresponding curve will yield a detectable signal. 
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FIG. 6: The flux of 77 monochromatic 7-rays as a function of m x , assuming (<tw) T7 ~ 2 x 10 30 
cm s , and for the SIS, the Moore et al., and the NFW profile. The clump masses used are 1O 2 M 
and 1O 8 M . Also shown as functions of _E 7 = m x , are the two dominant background contributions, 
the hadronic cosmic ray and the electronic cosmic ray induced ones. The dashed line represents 
the minimum flux, F m i n , required so that a 5-cr detection be achieved for A e ff = 10 8 cm 2 , E t h=50 
GeV, ag = 6', and 100 hrs of observation. For m x values corresponding to fluxes higher than F m i n , 
the clumps will be detectable in the 77 line. 
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FIG. 7: The minimum detectable (av)zj as a function of m x for the SIS, the Moore et al., and 
the NFW profile. The clump masses used are 1O 2 M0 and 1O 8 M0. The dots represent allowed 
SUSY models. We require a 5-cr detection level for A e ff = 10 8 cm 2 , £"^=50 GeV, ag = 6', energy 
resolution Re = 15%, and for 100 hrs of observation. For a specific clump mass and a specific 
density profile, and for the above instrument and observation characteristics, only the SUSY models 
that lie above the corresponding curve will yield a detectable signal. 
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FIG. 8: The flux of Z7 monochromatic 7-rays as a function of m x , assuming (av)zj — 3 x 10 30 
cm s , and for the SIS, the Moore et al., and the NFW profile. The clump masses used are 1O 2 M 
and 1O 8 M . Also shown, as functions of m x , are the two dominant background contributions, the 
hadronic cosmic ray and the electronic cosmic ray induced ones. The dashed line represents the 
minimum flux, F m i n , required so that 5-cr detection be achieved for A e ff = 10 8 cm 2 , E t h=50 GeV, 
<j0 = 6', and 100 hrs of observation. For m x values corresponding to fluxes higher than F m i n , the 
clumps will be detectable in the Z7 line. 
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